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Abstract 

We propose to observe QSO-galaxy strong lens systems to give a new constraint 
on the damping scale of the initial fluctuations. We find that the future observation 
of submilliarc scale astrometric shifts of the multiple lensed images of QSOs would 
find ~ 1O( 3 ~ 9 )M0 subhalos inside the macrolens halo. The superweakly interacting 
massive particles (superWIMPs) produced from a WIMP decay and the warm dark 
matter (WDM) particles that predict a comoving damping scale larger than ~ 2 kpc 
can be constrained if ~ 10 3 Mq subhalos are detected. 



1 Introduction 



There has been mounting evidence that most of the matter in the universe is not luminous 
but dark. Current observations such as the cosmic microwave background (CMB) suggest 
that the dark matter (DM) makes up about 25 % of the universe pp. The dark matter 
is usually assumed to be a collisionless cold component (non-relativistic at the time of 
freeze-out), called as the cold dark matter (CDM). Weakly-interacting massive particles 
(WIMPs), such as the lightest neutralino in the supersymmetric standard model (SUSY 
SM) j2], and the lightest Kaluza-Klein particle in the universal extra dimension (UED) 
[3], are the popular CDM candidates. The predicted thermal relic abundances and the 
large-scale structure (^ 1 Mpc) are in good agreement with the observed values. 

However, the recent high-resolution iV-body simulations on the CDM-based structure 
formation revealed various discrepancies on smaller scales (J$ 1 Mpc). The first one is so- 
called the "missing satellite problem" ^j: the iV-body simulations of the CDM particles 
predict significantly more virialized dark objects with mass M^z 1O 9 M (or subhalos) in 
the galaxy-sized halos with M ~ 1O 12 M than those observed around the Milky Way. The 
other one is called the "cusp problem" j^j: the CDM-based models also predict a cuspy 
profile for mass density distributions for the CDM halos [S] although the measurements of 
the rotation curves imply the presence of cores in the centers of the halo. 

Although these discrepancies may be circumvented by some baryonic processes jjj, it 
may be worthwhile to consider the other kinds of DM particles with different clustering 
properties. For instance, the superweakly interacting massive particles (superWIMPs) [8J 
or the warm dark matter (WDM) [§| particles can have large velocity dispersion at the 
epoch of radiation-matter equality. If the DM consists of the superWIMPs or the WDM 
particles, the number of less massive ^ 1O 9 M subhalos is significantly reduced, and the 
cusp formation is also suppressed because the primordial fluctuations at the small scales 
(^ 1 Mpc) are damped. 

In this note, we consider the possibilities of probing such subhalos with M ^ 1O 9 M 
via strong lensing in which the image separations are on submilliarcsecond scales, called 
as "submillilensing" . Recently it has been pointed out that the future submillilensing 
observations of multiply-imaged QSO-galaxy lens systems can directly probe the mass 
scale of subhalos in the parent galaxy halo ^U]. The submillilensing observation might 
resolve whether above problems originate from various baryonic contrivances or the nature 
of the DM particles. In section 2, we begin with the discussion of submillilensing and 
consider the possibility of detecting the small-mass subhalos via submillilensing in the 
next decade. In section 3, we study its implications to the superWIMP and the WDM 
scenarios. Section 4 is devoted to Conclusions and Discussion. 

2 Submillilensing 

In what follows, we discuss the possibilities of direct detection of subhalos with a mass of 
M ^ 1O 3 M via substructure lensing which is defined as lensing by M ^ 1O 9 M subhalos 
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that perturb a "simple" strong lensing by ~ 10 12 M Q parent galaxy halo. To date, about 
10 quadruply-imaged gravitational lenses with flux ratio "anomalies" have been detected 

IT2| IT3] . Here "anomaly" refers to an observed image flux ratio that does not agree with 
the ratio predicted by standard macrolens models with a smooth gravitational potential. 
From the radio and the mid-infrared observation, some of those lens systems are found to 
be consistent with the model in which the macrolens with a smooth potential is perturbed 
by subhalos [TTJ [TJJ EH EH Hi] . Unfortunately, the mass scale of the subhalos has not 
been yet determined well. 

First, we theoretically estimate the mass range of the subhalos that can perturb the 
fluxes of the multiple images produced by the parent galaxy halo with a mass of ~ 1O 12 M . 
In what follows, for simplicity, we assume that the low-mass subhalos are described by 
tidally-cut singular isothermal spheres (SISs) with a mass function dn/dm oc m~ 2+e , where 
< £ <C 1. At a distance r from the center, an SIS with a one-dimensional velocity 
dispersion ogis has a density profile p(r) oc Osis/ 7 "" 2 [TE] . If we further assume that the 
parent galaxy halo is also described by an SIS with a one- dimensional velocity dispersion 
a, then the tidal radius r t of the SIS subhalo with <tsis at a distance r from the center 
of the parent halo is approximately given by (Tsis/rt ~ cr/r, which yields the mass of the 
tidally-cut SIS as m oc r<7g Is . The effects of deviation from the assumption made here will 
be discussed later. 

The lensing cross section by a subhalo is proportional to the square of the Einstein 
angular radius 8e, and the radius for the SIS is proportional to <7g IS [16J. Using the 
equations for the mass function and the tidal radius, the substructure lensing cross section 
per logarithmic mass interval is given by 

dr oc 9 2 E dn oc 0g is dri ~ m 1//3+e <i(lnm) (2.1) 

where we have assumed that the distance r of the subhalo to the halo center is approx- 
imately equivalent to the Einstein radius of the parent halo r^. This approximation is 
verified because the substructure lensing cross section depends weakly on the distance to 
the center r as oc r _4//3 . The contribution from the distant halos in the line of sight can 
boost the cross section by a factor J7|. Thus, the contribution of massive subhalos is 
significant in the substructure lensing as long as the mass function satisfies e > —1/3. 
Assuming that the mass function with < e < 0.2, observed for the massive scales 
M > 1O 9 M [18J, is applicable to the mass scale 1O 3 M < M ^ 1O 9 M , the lensing proba- 
bility by an intermediate-mass subhalo ~ 1O 3 M is reduced by ~ 10~( 3 ~ 4 ) in comparison 
with the probability by a massive subhalo ~ 1O 9 M . 

If the subhalos have a different mass profile such as the NFW profile jE], the contri- 
bution of less massive subhalos can be further reduced. Because the ratio between the 
Einstein radius of the perturber and the tidal radius at a fixed distance r decreases as the 
subhalo mass decreases, the lensing cross section for a less massive subhalo is susceptible 
to the inner less cuspy profile. Thus, we expect that the contribution from M ^ 1O 3 M 
subhalos is negligible in alternating the amplitude of the flux of multiple images. 
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Next, we explore the possibility of measuring the mass scale of these subhalos from 
astrometric shifts of the multiply lensed images. For SISs, the order of the magnitude of 
the astrometric shifts is typically given by the size of the angular Einstein radius 9e [19] - 
For an SIS with one-dimensional velocity dispersion asis at ~Gpc, 9g is approximately 
given by 



Thus, observation with angular resolution of submilliarcsecond scales (~ O.Olmas), which 
will be achieved in the next generation satellite VLBI mission such as the VSOP2 [20], can 
reveal subhalos with one-dimensional velocity dispersion asis > 0.6 km s _1 . It corresponds 
to M^1O 3 M at the distance equal to the Einstein radius r = te of the macrolens, 
assuming that the velocity dispersion of the parent halo is a ~ 200 km s _1 . 

From the astrometric shifts of the multiple extended images perturbed locally by a 
subhalo with respect to an unperturbed macrolensed image, we can break the degeneracy 
between the subhalo mass and the distance in the line of sight to the images if resolved at 
scale of an Einstein radius of the perturber ^TJEJ]. This is of great importance because 
otherwise we cannot determine whether the flux ratio anomaly is caused by more massive 
intergalactic halos in the line of sight or by less massive subhalos within the macrolens 
halo. Even if the density profile of the perturber is shallower than an SIS, we can make 
a distinction between models with different density profiles from astrometric shifts of 
the surface brightness profile within the source [IHj. A direct detection of less massive 
10 3 M Q ^5 M ^ 1O 9 M subhalos within the parent halo will give a stringent constraint on 
the superWIMP and the WDM scenarios, which will be discussed in the next section. 

3 Implications to SuperWIMP and WDM Scenarios 

There exist many well motivated models for superWIMPs and WDM from particle physics. 
The natural candidates for the superWIMPs are gravitino and axino, that are superpart- 
ners of graviton and axion, respectively [H|. The right-handed sneutrino is also the can- 
didate when neutrino masses are Dirac-type [21]. Others are Kaluza-Klein graviton and 
axion states in the UED jSJ. As for the WDM, a light gravitino and sterile neutrinos 
have been discussed as such candidates. Thus, it is interesting to study the feasibility of 
probing or constraining those models in the near future experiments. 

In the following, we refer "superWIMPs" as the particles whose interactions are weaker 
than the weak interaction, such as the gravitino which couples to other particles only 
gravitationally. SuperWIMPs are also assumed to be produced by the decays of heavier 
particles whose interactions are weak {e.g. WIMPs)* 1 . 

In the CDM-based structure formation models, the structure of the universe forms 
in a hierarchical manner. Protohalos, which are the first virialized objects, appear first 

^Usually conventional WDM particles do not satisfy both of these properties. In this sense, super- 
WIMPs can be distinguished from conventional WDM particles. 




(2.2) 
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after the mass density fluctuation becomes nonlinear, and larger objects form successively 
via their merger. In the ordinary WIMP models, the comoving damping scale of the 
power spectrum is typically (0.01 — 10)pc, depending on the Hubble radius at the kinetic 
decoupling temperature |22j, and the protohalo mass is (10~ 12 — 1O~ 4 )M [23J. 

In the superWIMP or WDM scenarios, the primordial fluctuations on even larger scales 
can be damped because they have large velocity dispersion at the decoupling. Therefore, 
the protohalos become massive, and the total amount of subhalo mass inside the parent 
halo can be significantly reduced. 

When the comoving damping scale of the DM power spectrum is R cn t, the protohalo 
mass is roughly given by 



Here, the present matter density f2 m ,o = 0.24 and the Hubble constant h = 0.73 are 
assumed [T]. In the previous section, we showed that the sensitivity of direct detection of 
subhalos inside the parent halo with submillilensing may reach intermediate mass ~ 1O 3 M 
scales. Thus, when R cut ^ 2 kpc, the future submillilensing experiments may directly 
detect such protohalos. 

First, we discuss the damping scale in the superWIMP scenario and compare it with 
the sensitivities of future submillilensing. In the scenario, some WIMPs freeze out from the 
thermal equilibrium as the usual WIMP DM models, and superWIMPs are non-thermally 
produced from the WIMP decay, since the superWIMPs interact superweakly with the 
thermal bath. The scenario retains the property of the ordinary WIMPs that the observed 
relic density is naturally achieved. Furthermore, they are produced by decay of some long- 
lived massive particle X. Since they can have large velocities at the production epoch, 
the free-streaming damps the small-scale inhomogeneities. The superWIMP scenario is a 
possible explanation for the small-scale structure |2*4^ I25j. 

The comoving damping scale is typically given by the free-streaming length of the 
superWIMP, R {s , at the matter-radiation equality t eq . When the long-lived particle has 
lifetime Tx and the superWIMP is produced with three-momentum normalized by its mass 
u(= p/m), the comoving free-streaming scale Rf s is given by 



where w eq = (a{jx) / 'a(t cq ))u and a(t) is the scale factor as a function of t. 

When the superWIMPs are produced from decay of electrically-charged particles* 2 , 
the small-scale power may be further suppressed 28;. The charged particles are coupled 

# 2 It is recently pointed out in Ref. [2H] that long-lived charged particles (tx ^ 10 3 sec) lead to overpro- 
duction of 6 Li in the Big Bang Nucleosynthesis due to the catalytic enhancement of the production. On 
the other hand, it is argued in Ref. that there are a lot of ambiguities in the derivation. 




(3.3) 
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to the photon-baryon fluid which oscillates on subhorizon scale. If the scale in question 
enters the horizon before the decay, the density fluctuation of such DMs cannot grow but 
oscillates, which gives a rise to the suppression on small scales. The damping scale is 
typically given by the Hubble radius at the decay time, H^ijx), 



In Fig. 1 we show the damping scales R cut of the charged and neutral X 
functions of Tx and u. When X is neutral, the damping scale is determined by the 
free streaming scale Rf s . Longer lifetime tx implies a larger free streaming scale, since 
the velocity dispersion at the radiation-matter equality time, which is proportional to 
a(rx)/a(t e q), becomes larger. When X is charged, the damping scale R cut corresponds to 
the larger one of _Rf s and R c ^. For u <C 1, R cut is given by i? c h, since i?f s is ~ u x R ch . 
The hatched region with damping scales larger than ~ 1 Mpc is constrained from Lyman 
alpha clouds The future submillilensing experiments may cover regions above the 

bold lines, which correspond to R cut ^ 2 kpc. When the X lifetime is longer than ~ 1 sec 
and the super WIMPs are produced with relativistic momentum (u ^ 1), the damping scale 
is larger than ~ 1 kpc, which may be constrained by future observation if the subhalos 
with M ~ 1O 3 M were discovered. When X is charged, the region with t x > 400 sec may 
be also covered even in the small u cases. 

Following Ref. [21], we indicated the parameter region which is suitable to solve the 
discrepancies in small-scale structure in Fig. 1. The region with R cnt ~ (0.4 — 1.0) 
Mpc is suitable to solve the "missing satellite problem", while the gray region with 
1 ~ n _3 (rx/10 6 sec) 3/ ' 2 ^ 4 is favored to solve the "cusp problem", which requires a large 
DM velocity dispersion [30J. As we can see in Fig. 1, the model parameters corresponding 
to these regions can be well constrained by the future submillilensing observation. 

For the illustrative purpose, we consider a superWIMP gravitino model, in which the 
gravitinos are produced by slepton or sneutrino decay. Sleptons or sneutrinos are assumed 
to be the lightest SUSY particles in the SUSY SM* 3 . The slepton and sneutrino lifetimes 
are given by (3^] 



where = 2.4 x 10 18 GeV, m is the gravitino mass, and mx is the mass of the parent 
WIMP X. In Fig. 2, R^ s (solid lines) is shown as a function of m and the mass difference 
between the gravitino and the parent particle (Am = mx — rn). Dashed lines indicate 
the damping scale for which R c ^ > Rf s . The gray region which corresponds to explain the 

# 3 Thc lightest ncutralino is also one of the candidates for the lightest SUSY particle in the SUSY SM. 
However, the decay produces hadronic shower to spoil the BBN if it is not photino-like. The long-lived 
slepton and sneutrino are not strongly constrained by the BBN, since their hadronic branching ratios are 
small. See Ref. |3J for the constraints on the superWIMP gravitino model. 
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"cusp problem" [21] m ay be excluded when the future submillilensing experiments find 
subhalos with mass smaller than ~ 1 x 10^ 7_8 ^M Q . It is known that the lifetime and the 
mass of the long-lived particles whose decay produces the superWIMPs, are constrained 
by the Big Bang Nucleosynthesis (BBN) and by the CMB Planckian spectrum, depending 
on the decay channels. When the hadronic modes are dominant in the decay, the BBN 
constrains the lifetime to be shorter than ~ 1 sec [33]. Even if the hadronic shower is 
suppressed, the electromagnetic energy injection from the long-lived particle decay to the 
thermal bath is also constrained from the BBN while the constraint is weaker. In addition 
to it, the CMB Planckian spectrum also gives a constraint when the lifetime is longer than 
~ 10 6 sec [3J. 

The future submillilensing experiments are complementary to those constraints, since 
the damping scale R cut is independent of the decay channels. For example, when the 
superWIMP gravitinos are produced from sneutrino decay, the energy injection to thermal 
bath is very tiny so that the constraints from the BBN and the CMB Planckian spectrum 
are very weak. Even in such a case, the submillilensing will still constrain the model. 

Next, we discuss the damping scale in the WDM models. Among many WDM candi- 
dates discussed so far, the light gravitino with mass of order from 10 -6 eV up to 1 keV, 
which is likely to be the lightest SUSY particle in gauge-mediated models of supersymme- 
try breaking, is one of the well-motivated models from the viewpoint of particle physics 
[3U EE]- Such gravitinos are in thermal equilibrium at early times but decouples when the 
degrees of freedom g*{Tr>) is 0(100) where To is the decoupling temperature. Because the 
decoupling temperature of such species is higher than that of (active) massive neutrinos 
which play roles of hot dark matter, their velocity dispersion is not so large compared to 
that of massive neutrinos but non-negligible at the time of structure formation. Thus, 
they can act as the WDM. The comoving damping scale for the free-streaming for such 
gravitinos or any other thermal relic can be written as 

Rk „ o.84Mpc (»m (^T) (Mil) , (3.7) 
1 V 10-75 ) V"»wdm/ V 315 ) ' 

where (p/T) is the mean momentum over the temperature. For thermally decoupled 
species, this factor gives almost unity, i.e., (p/T) /3.15 ~ 1. The requirement from Lyman 
alpha clouds, i?f s ^ IMpc, implies mwDM ~ IkeV. On the other hand, the energy density 
of WDM can be written as O dm /j 2 = (mwDM/94eV) (10.75/g* (T D ) ). Assuming Q DM h 2 ~ 
0.10, the mass of WDM should be ijiwdm ~ 0.1 keV even for g*(To) ~ 100. We need 
to introduce more extra degree of freedom around Tp as g*(To) ~ O(10 3 ). When the 
constraint on Rf s is improved to be ^5 lkpc, the lower bound for the mass can reach 
wwdm ~ 1 MeV for g*(Tr>) ~ 100. Thus, the WDM scenario may face further difficulties 
if future submillilensing experiments would find small-mass subhalos. 

Another well- motivated candidate for WDM is the sterile neutrinos [31] • Because they 
directly couple to the active neutrinos alone, they can be produced via neutrino oscillation. 
Although the evaluation of their energy density requires a numerical integration of the 
Boltzmann equation, some useful fitting formulae are available. The present energy density 
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of sterile neutrinos can be written as [37] 



n Us h 2 ~ 0.3 



( 



sin 2 29 



) 



lOOkeV 



m 



s 



) 



2 



(3.8) 



where 6 is the mixing angle between the active neutrinos and sterile neutrinos and m s is 
the mass of sterile neutrinos. The temperature at the time when the production is most 
efficient is jSZl EHJ 



The sterile neutrinos can damp the small-scale inhomogeneities via the free streaming 
in the same manner as the thermally decoupled WDM particles do. The free-streaming 
scale for a sterile neutrino WDM can also be obtained using Eq. ()3.7|) with different values 
for (p/T) from that of the thermally decoupled ones. Because the sterile neutrinos are 
not in the thermal bath at early times, their distribution function deviates from that of 
a thermal one and the above factor can be (p/T)/3.15 ~ 0.9 for the standard production 
mechanism jSH]. Although there are some differences between the thermally decoupled 
WDM and the sterile neutrino WDM, they give the same predictions for the damping of 
matter power spectrum by identifying their masses as [""""H I"""!] 



where m t h er mai is the mass of thermally decoupled relics such as a light gravitino, which is 
denoted as m WDM in Eq. (|3.7|) . Because the shape of matter power spectrum is determined 
by the ratio m/T and the density parameter f2rjM^ 2 - Accordingly, the constraint on the 
sterile neutrino mass bound would be different from that on the thermally decoupled WDM 
mass bound* 4 . Using the above formula and assuming the damping scale as R cut ~ 1 kpc 
which can be reached by future submillilensing experiments, we can expect that the mass 
of sterile neutrinos can be constrained to be m s ^ 40 MeV, which will be in conflict even 
with the current constraint m s ^ 10 keV [3*9*1 E21 ESI- Thus, we can obtain much more 
insight on the robustness of these scenarios from the future submillilensing experiments, 
as well as for the superWIMP scenarios. 

Here some comments on the mixed dark matter scenario are in order. It is possible 
that, for example, gravitinos are produced not only from the decay of the next-lightest 
supersymmetric particles (NLSP) but also from thermal plasma. In this case, the present 
DM is composed of CDM and superWIMP. In such mixed DM scenarios, the damping of 
the small-scale structure is less significant in comparison with the models in which the 
superWIMPs make up all of the DM. To what extent the amplitude at small scales can 
be reduced depends on the ratio of the energy density of CDM and superWIMP DM. 
Detailed discussion on this issue is beyond the scope of this letter. Some discussions on 

# 4 The properties of the sterile neutrinos can also be constrained from the measurement of the X-ray 
flux, since they can contribute to the X-ray flux due to radiative decay [52]. See also Refs. [3*2 [""""]• 




(3.9) 




(3.10) 
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the matter power spectrum in such mixed models can be found in Refs. jH] and 03] . |Hj 
analyzed models in which the superWIMPs are produced from the charged particle decay 
while 01] considered models in which the superWIMPs are produced from the neutral 
particle decay. 

4 Conclusions and Discussion 

We have shown that future observation of multiply-imaged QSO-galaxy lens systems with 
a high angular resolution ~ 0.01 mas will prove the small-scale clustering properties of the 
DM halos down to ~ 1O 3 M . The presence of ~ 1O 3 M subhalos implies the comoving 
damping scale of the primordial fluctuations ~ 2 kpc assuming that the observed subhalos 
retain the original mass during the merger process. The superWIMP and the warm DM 
scenarios that predict a larger damping scale (0.002—1) Mpc in comparison with (0.01 — 10) 
pc in the ordinary WIMP scenarios would be strongly constrained if the presence of such 
subhalos were proved. 

In the superWIMP scenario, the superWIMP DM is produced from the long-lived 
particle X decay, and free-streaming of the superWIMPs damps the small-scale inhomo- 
geneities. When the lifetime of X is longer than ~ lsec, the damping scale is larger than 
~ 1 kpc unless the superWIMP and X masses are degenerate and the superWIMPs are 
non-relativistic at the production epoch. In addition, when X is a charged particle, it is 
coupled to the oscillating photon-baryon fluid before the decay. Therefore, the small-scale 
inhomogeneities inside the sound horizon cannot grow. The damping scale becomes larger 
than ~ 1 kpcx (rx/100sec) 1,/2 in the case of charged X. One of the natural superWIMP 
candidates is the gravitino produced from the slepton or sneutrino. In this case, the 
damping scale is typically larger than 10 2 kpc. The future submillilensing experiments, 
which cover the DM subhalos with mass ^ 1O 3 M , are important tests for probing such 
superWIMP scenarios. For the WDM scenarios, such as the light gravitino and the sterile 
neutrinos, the WDM mass would be further constrained. 

The survivability of the protohalos during the merger process is under debate now. 
In the ordinary WIMP scenarios, it is claimed that most of the earth-mass protohalos 
are stable against the tidal stripping [46J. It is also discussed whether those halos are 
disrupted by interaction with stars 07] . On the one hand, subhalos that cross the galactic 
disk nearly perpendicularly or that fall off to the center of the parent galaxy are strongly 
disrupted by the tidal force, leading to a significant decrease in the total mass. On the 
other hand, subhalos that orbit on the plane nearly parallel to the disk or that reside in 
the low density region survive more-or-less intact. 

For simplicity, we have assumed that the observed mass scale of the subhalo is equiva- 
lent to the lowerbound on the protohalo mass scale. In practice, however, the observable 
mass using gravitational lensing is limited to the one within a certain radius centered at 
the line of sight. As a result, there remains a certain ambiguity in estimating the mass 
scale of the observed subhalo. Observation of astrometric shifts of lensed QSO images 
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with a substructure in the surface brightness may help to reconstruct the subhalo mass 
density profile, thereby reducing the ambiguity [TUj . 

It has been argued that the superWIMP and WDM scenarios can resolve the small-scale 
^5 1 Mpc discrepancies, such as the "missing satellite problem" and the "cusp problem" 
if the damping scale is as large as (0.4 — 1.0) Mpc. Therefore, future submillilensing 
experiments will shed a new light on these small-scale structure problems once the above 
ambiguities are removed. 
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Figure 1: Damping scales of the DM power spectrum as functions of Tx and u in the 
superWIMP scenarios. The upper figure is for the neutral X case, and the lower one 
is for the charged X one. The (hatched) region with damping scale larger than 1 Mpc 
is constrained from Lyman alpha clouds. The gray region is favored to solve the cusp 
problem. The future submillilensing observations may cover the region above the bold 
lines which correspond to R cut ^ 2 kpc. 
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Figure 2: Damping scales for the gravitino superWIMPs produced by slepton or sneutrino 
decay. Rf s (solid lines), (dashed lines) are plotted as functions of m and Am = mx — m, 
respectively i? c h is shown only for the case _R c h > i?f s . The gray region is favored to solve 
the cusp problem. 
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